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We study the phenomenon of parton-hadron duality in both polarised and unpolarised electron 
proton scattering using the HERMES and the Jefferson Lab data, respectively. In both cases we 
extend a systematic perturbative QCD based analysis to the integrals of the structure functions in 
the resonance region. After subtracting target mass corrections and large x resummation effects, 
we extract the remaining power corrections up to order 1/Q^. We find a sizeable suppression of 
these terms with respect to analyses using deep inelastic scattering data. The suppression appears 
consistently in both polarised and unpolarised data, except for the low polarised data, where 
a large negative higher twist contribution remains. Possible scenarios generating this behavior are 
discussed. 

PACS numbers: 13.88.+e, 13.60.Hb 



I. INTRODUCTION 

The structure of hadrons and their interactions can 
be described within two different but complementary ap- 
proaches based on either partonic or hadronic degrees 
of freedom. The first one is expected to be valid at 
high energy, while the second one is applicable at low 
energy where the effects of confinement become large. 
In some specific cases where a description in terms of 
non-partonic degrees of freedom seems more natural, the 
quark-gluon description can be also successfully used. 
This observation is called parton-hadron duality. It was 
introduced for deep inelastic scattering (DIS) by Bloom 
and Gilman Q who reported an equivalence between the 
smooth X dependence of the inclusive structure function 
at large and the average over W'^ of the nucleon reso- 
nances {x = /2Mv, is the four momentum transfer 
squared, M is the nucleon mass, v is the energy trans- 
fer, and W"^ = Q'^il/x - 1) + is the final state in- 
variant mass). One refers to global duality if the aver- 
age, defined e.g. as the integral of the structure func- 
tions, is taken over the over the whole resonance region 
1 < W'^ < 4 GeV^. If, however, the averaging is per- 
formed over smaller ranges, extending e.g. over sin- 
gle resonances, one can analyze the onset of local duality. 

More generally, the concept of duality is often assumed 
in QCD-based interpretations of most hard scattering ex- 
periments, such as DIS, e^e~ annihilation into hadrons, 
and hadron-hadron collisions. Its usage appears when- 
ever hadronic observables (mostly averaged over a given 
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energy range) are replaced by calculable partonic ones 
with little more going into the hadronic formation phase 
of each process - from partons to hadrons or vice versa. 
In a phenomenological context, duality studies are aimed 
at establishing to what extent a partonic description of 
the hard scattering process can determine the structure 
of the final state. In fact, as prescribed by the factoriza- 
tion property of QCD, we visualize hard scattering pro- 
cesses happening in two stages, one dominated by short 
times and distances and involving only parton jets, fol- 
lowed by hadron formation at a much larger scale. Du- 
ality is intrinsic to the factorization property. Violations 
of duality might signal violations of factorization in that, 
for instance, the probe-parton interaction might occur 
at larger time scales than required in order to exclude 
parton (re)interactions. 

With the advent of both more detailed studies of soft 
scales and confinement , and higher precision measure- 
ments covering a wide range of reactions, it is now be- 
coming possible to investigate the role of duality in QCD 
as a subject per se. For example, recent studies of local 
parton-hadron duality and its violations in semi-leptonic 
decays, and t decays illustrate how the possible impact 
of these experiments on the extraction of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements, depends 
on the size of violations of local duality 01 ■ A practi- 
cal necessity to address duality quantitatively exists also 
for inclusive ep scattering where most of the currently 
available large x data lie in the resonance region. In fact, 
for X > 0.5 and ^ 5 GeV^ - a typical starting value 
for perturbative evolution - < 5 GeV^. Therefore, 
the behavior of the nucleon structure functions in the 
resonance region needs to be addressed in detail in order 
to be able to discuss theoretical predictions in the limit 
X —* I. 

The first QCD-bascd studies of Bloom and Gilman du- 
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ality reinterpreted the "averaging" procedure in terms 
of Mellin moments of the structure function. The mo- 
ments taken in the low and in the DIS regime, re- 
spectively, were shown to be equivalent to one another 
within the given range and precision of the data, mod- 
ulo perturbative corrections and relatively small power 
corrections Q. It was conjectured that duality resulted 
from a cancellation of higher order terms in the twist 
expansion that would otherwise be expected to domi- 
nate the cross section at x — > 1. This view has been 
adopted since, particularly in the more recent studies in 
Ref. In Ref. a new analysis was performed, using 
the recent inclusive unpolarised electron-nucleon scatter- 
ing data on hydrogen and deuterium targets from Jef- 
ferson Lab It was shown in particular that, because 
of the increased precision of the data, one is now able 
to unravel different sources of scaling violations affect- 
ing the structure functions, namely Target Mass Correc- 
tions (TMC), Large x Resummation effects (LxR), and 
dynamical Higher Twists (HTs) , in addition to the stan- 
dard Next-to-Leading-Order (NLO) perturbative evolu- 
tion. As a result, contrarily to what originally deduced 
in e.g. Ref.Q, a more pronounced role of the HT terms 
is obtained, pointing at the fact that duality, defined on 
the basis of a dominance of single parton scattering, could 
indeed be broken. 

In contrast to the extensive study of duality for the un- 
polarised, i.e. spin averaged, photo-absorption cross sec- 
tion, the validity of duality has not been investigated un- 
til very recently for the spin structure function gi , which 
is proportional to the s^in- dependent photo-absorption 
cross section. Evidence of duality for the spin asymmetry 
Ai was reported in Ref. A phenomenological study 
addressing parton-duality was performed in Q using the 
low data from Ref. [13 ■ Studies of duality for gi are of 
particular interest because they might help understand- 
ing the transition from the large regime described 
by pQCD, and the ^ limit, where the Gerasimov- 
Drell-Hearn sum rule is expected to apply They 
may also lead to a complementary method to study the 
spin structure of the nucleon at large x, which is diffi- 
cult to measure in the DIS region with high statistics. In 
particular, they might provide additional information on 
the transition from single parton scattering, to the dom- 
inance of processes where several partons are involved 
|l2l |. In this respect, it is important to perform an anal- 
ysis aimed at disentangling the different contributions to 
the dependence of gi in the resonance region. The 
aim of this paper is to carry out such an analysis by in- 
vestigating quantitatively the onset of duality and its vi- 
olations both for the unpolarised and polarised structure 
functions. 

In Section II we define the concept of duality and 
we illustrate the role of different kinematical regions; 
in Section III we present our analysis and we describe 
our results. In particular, we compare the data with 
perturbative-QCD predictions and we discuss in detail 
the contribution of different types of corrections in both 



the unpolarised and polarised case. Finally, in Section 
IV we draw our conclusions. 



II. DEFINITIONS AND KINEMATICS 

Parton-hadron duality in DIS was first observed more 
than three decades ago. Since then it has been necessary 
to give new definitions of the quantities involved which 
can be described within QCD-based approaches. In what 
follows we list all such definitions. 



A. Kinematical variables 

Besides the scaling variable x, other variables have 
been used in the literature to study duality. A num- 
ber of parameterisations based on these variables have 
been proposed that reproduce in an effective way some 
of the corrections to the perturbative QCD calculations 
that we study in this paper. The most extensively used 
variables are: x' = \/lo' , where oj' = 1/x -\- /Q^. x' 
was originally introduced by Bloom and Oilman in or- 
der to obtain a better agreement between DIS and the 
resonance region; ^ = 2a;/(l + [l + ^x^W^ IQ'^fl'^) pl| . 
originally introduced to take into account the target mass 
effects; x^ = Q"^ \- B/{Q'^ + - M"^ + AYA and B 
being fitted parameters, used in Refs. 0, 0. These 
additional variables include a dependence that phe- 
nomenologically absorbs some of the scaling violations 
that are important at low Q^. In FigQlwe compare their 
behavior vs. x for different values of . From the figure 
one can see that by calculating F2 in ^ and x' , one effec- 
tively "rescales" the structure function to lower values of 
X, in a dependent way, namely the rescaling is larger 
at lower Q^. 

In this paper we present results in terms of x and 
and we illustrate the contributions of scaling violations 
of different origin on a case by case basis. 

B. Unpolarised structure function. 

The inclusive DIS cross section of unpolarised electrons 
off an unpolarised proton is written in terms of the two 
structure functions F2 and Fi, 

d^a _ 47rQ^ r/ ^ {Mxy)^ \^ ^ v^xf] (1) 
dxdy Q^xy [\ J J ' 

where y = h'/ E, E being the initial electron energy. The 
structure functions are related by the equation: 

Fi =F2(l+7')/(2x(l + i?)), (2) 

where 7^ = A^M^x^ jQ'^ ^ and R is ratio of the longitudinal 
to transverse virtual photo- absorbtion cross sections. 

In QCD, Fi is expanded in series of inverse powers 
of . obtained by ordering the matrix elements in the 
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DIS process by increasing twist r. which is equal to their 
dimension minus spin: 



is unity. However, one can extend this definition to the 
ratio: 



+ 0(1/Q4) (3) 



R 



HT 



JHT ■ 



(9) 



The first term is the leading twist (LT), with r = 2. The 
terms of order t > 4, in Eq.Q are the higher 

order terms, generally referred to as higher twists (HTs). 
Additional corrections to the LT part due to the finite 
mass of the initial nucleon - the target mass corrections 
(TMC) - are included directly in i^2^^. For larger 
than « 1 GeV^, TMC are taken into account through 
the following expansion [T^ : 



The rational for this definition is that even when in- 
cluding the first few terms of the twist expansion, one is 
embracing a partonic description of the proton. 

Perturbative QCD analyses use the Mellin moments of 
the structure function, which allow for a direct compar- 
ison with theoretical predictions. These are defined as: 



LT(TMC) 



(4) 



Jo 



(10) 



and by 



where F.^'^^ is the structure function in the absence of 
TMC. Since TMC should in principle be applied also to 
the HT, we disregard terms of 0(1/(3") 0. H, then, 
represents the "genuine" HT correction that involves in- 
teractions between the struck parton and the spectators 
or, formally, multi-parton correlation functions. 

Parton-hadron duality in DIS is studied by considering 
integrals of the structure function defined as: 



6(n-l) /f 



(7i + 2)(n + 3) 



n{n-\) (i 



^-1 ' (12) 
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Fl''\x,Q') dx 



(5) 



where F^®** is evaluated using the experimental data 
in the resonance region. For each value: x,nin = 

Wmin and Wmax delimit the resonance region. The same 
expression is then calculated in the same range of x and 
for the same value of Q^, using parameterisations of F2 
that reproduce the DIS behaviour of the data at large 
Q^. These parameterisations are very well constrained 
in the region of interest {x > 0.3) although they do not 
correspond directly to measured data. Here in fact F2 is 
dominated by the valence contribution. On the contrary, 
by using the same procedure at low x where the singlet 
and gluon distributions govern F2, one would find much 
larger uncertainties in the initial low parameterisa- 
tions because of their strong correlation with the value 
of as- We present two forms for the DIS integrals: 



and 



Fi^{x,Q^)dx, 



F2'^{x,Q') + ^^^^ 



(6) 



dx. (7) 



Duality is attained strictly only if the ratio: 

^rcs 

Jlt ' 



TjLT 
unpol 



(8) 



{n + 2){n + 3) \x 

which takes into account TMC However, one needs 
in this case experimental values of the structure func- 
tion in kinematics outside the resonance region. This 
procedure renders the comparison between theory and 
experiment less straightforward. The difference between 
the Mellin moments fEas llOll^ and the integrals over 
the resonance region - /„ = a;"~^F^"^(a;, Q^) dx - 

is shown in Fig|21 The drop of the quantities /„ with 
respect to M„ at larger values of Q^, is due to the pQCD 
evolution of F2, that moves strength to lower values of x, 
outside the range [smax, a^min]- In our approach we use 
the integrals defined in Eas lfcil7l As it can be understood 
also from the trend shown in Fig|2 these effectively de- 
scribe duality also as a function of the average value of x 
in each interval [a;min(g^), 2:max(Q^)]- This corresponds 
to: (a;) = x{W'^ = 2.5 GeV^). 

It is also possible to consider a third approach, namely 
a point by point comparison of F2 both in the DIS and 
in the resonance region. The latter can in fact be fitted 
to a smooth curve tracing the resonances with a very 
high accuracy given by the increased precision of the new 
Jefferson Lab measurements. They are then compared 
directly to DIS parameterisations F2 at the same x and 
values. This procedure |^ is much less sensitive to 
the elastic contribution. All of the approaches described 
in this Section are necessary to interpret quantitatively 
the dependence of duality and of its violation. 



C. Polarised structure function. 

The spin-dependent part of the polarised deep inelastic 
cross section is given by: 
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dxdy 27r2Q' 
1-^-^7 I. 91 (a;,': 



(13) 



where _ffe and _ffjv ai'c the polarisations of the incident 
electron and of the nucleon of the target, respectively. 

Because of the mixing of gi and 52, a precise determi- 
nation of gi from a longitudinally polarised target alone 
is not possible. The experimentally measured cross sec- 
tion asymmetries are the longitudinal A\\ and the trans- 
verse ones, formed from combining data with oppo- 
site beam helicity: 



^,1 



rTT 



cri~* -I- crT" 



(14) 



The polarised structure functions are determined from 
these asymmetries: 



gi{x,Q'^) = 



Au + tan6'/2 ■ 



52(2;, Q^) = 



2d' 



E + E' cos e 
E' sin 61 



A, - All 



(15) 



where E' is the scattered Icpton energy. 6 is the scattering 
angle, d' = [(1 - e) (2 - y)] / [y (1 + eR{x[ Q^)] ^ with e being 
the degree of transverse polarization of the virtual photon 
and defined as e^^ = 1 + 2(1 + 7"^) tan2(6'/2). 

The virtual photon- absorption asymmetries Ai and A2 
are related to the measured asymmetries by: 



A\ 
A, 



D{A^ 
^ d{A2 - 



(16) 



where D is the depolarization factor D = y(2 — y){l + 
7V2)/[y2(l +7^)(l-2m ^/Q2)+2(l-y-7V/4)(l+i?)], 
and d = i:)v'2e/(l + e), ?? = 27(1 - y)/{2 - y) and C = 
77(1 -I- e)/2e are kinematic factors. 

From the measured asymmetries A\\ and A^, the vir- 
tual photon asymmetries can be related to the photon 
absorption cross section of the nucleon for a given x and 



the longitudinal asymmetry is measured, it is necessary 
to make an assumption for the asymmetry A2. 

From the measured asymmetry it is possible to evalu- 
ate the polarised structure function 91 ; in fact neglecting 
the term 7'^g2, Eg. ((171) reduces to: 



gi(x,Q2) « Ai(x)^^i(x,Q2). 



(18) 



This approximation is adequate for describing the data 
within their given accuracy. 

As for the unpolarised case, the twist expansion reads: 



g^{x,Q')=g^'^{x,Q') + 



H{x,Q^) 



0(1/Q4), (19) 



where, using Eas. H2l4ll8|l : 

5^(2^,(3') = Fi^ix,Q')Ar''{x) 

x{l + j^)/{2x{l + R'^P)), (20) 

where by the superscript "exp" we emphasize that we 
have used the experimental values for the quantities un- 
der consideration. 

We study parton-hadron duality by defining the inte- 
grals: 



-prcs 



gr{x,Q')dx, 



(21) 



where g™'' is obtained from the data in the resonance 
region, and 



J- 1 — 



g^^{x,Q'^)dx (22a) 



g^''{x,Q') + ^^^^)dx (22b) 



The ratios are given by: 



pres 
LT _ ^ 1 

pol 

^ 1 



nHT _ ^ 1 



pol ^HT ■ 

^ 1 



(23) 



As for the unpolarised case, duality is verified if either 
ratio is unity. 



Ai 



A. 



_ cri/2 - 0-3/2 
0'l/2 + 0-3/2 

CTl/2 + £^3/2 



C^LT _ 



91 



1^92 



Fi 

7(31+52) 



Fi 



(17) 



Here cri/2 and (T3/2 are the virtual photo-absorption 
cross section when the projection of the angular mo- 
mentum of the photon-nucleon system along the inci- 
dent photon direction is 1/2 or 3/2 respectively, ult is 
the interference term between the transverse and longi- 
tudinal photon-nucleon amplitudes, respectively given by 
CTT = (a'i/2 +0-3/2 )/2 and cttt = (0-1/2 -o-3/2)/2. If only 



III. ANALYSIS AND INTERPRETATION OF 
DATA 

In this Section we present a quantitative analysis of 
the Q" dependence of parton-hadron duality in both po- 
larised and unpolarised ep scattering. We take into ac- 
count all current data in the resonance region, 1 < < 
4 GeV^ . For the unpolarised case we used data obtained 
at Jefferson Lab in the range 0.3 < g2 < 5 GeV^ 0, and 
data from SLAG (^3| and references therein) for > 4 
GeV^. For the polarised case there are only few experi- 
mental data in the resonance region. One set is part of 
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the E143 data and it corresponds to — 0.5 and 
1.2 GeV^. Another set is the one from HERMES 
in the range 1.2 < < 12 GeV^. 

In the polarised case the dependence originates 
from the structure function Fi and from the ratio R.. 
In the evaluation of the denominators of Eas. l|23() . we 
used the SLAG global analysis |23 parameterisation for 
R. and, for Ai, a power law fit to the world DIS data 
at X >0.3: Ai = as already shown in Ref. This 
parameterisation of Ai is constrained to 1 at x=l and 
it does not depend on Q^, as indicated by experimental 
data in this range [2ll |. 



A. Comparison with pQCD 

The unpolarised structure function F2 was first eval- 
uated from dynamical parameterisations. coming from 
the Parton Distribution Functions (PDFs): MRST [2||, 
CTEQ m, GRV94 Q and GRV98 [H. The last two 
parameterisations have been evaluated at Leading Order 
(LO) and at Next to Leading Order (NLO). All of them 
are pure DIS parameterisations and they were extended 
to the measured x and ranges by pQGD evolution. 
The evolution of the polarised parton densities is gov- 
erned by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) 01 equations. The results are shown in Fig|31 
in the top panel for the unpolarised data and in the bot- 
tom panel for the polarised data. The uncertainty due 
to the use of different LO and NLO parameterisations 
is represented with a band; the other band represents 
the experimental uncertainty, calculated as the sum in 
quadrature of the statistical and systematic uncertain- 
ties of the data in the resonance region. 

Because the structure functions arc dominated by large 
X kinematics, DGLAP evolution proceeds only through 
Non-Singlet (NS) distributions. This explains why there 
is very little uncertainty in the extrapolation of the initial 
pQGD distribution to the low values of considered 
and the small difference between LO and NLO evolution. 

Parton-hadron duality is not fulfilled by using solely 
the PDFs up to NLO in both the unpolarised and po- 
larised structure functions F2 and gi . However it is pos- 
sible to see a different behavior between i?unpoi and i?poi. 
In the unpolarised case the ratio is increasing with Q^, 
but for the polarised case the situation is different: while 
at low the ratio is significantly below unity and shows 
a strong increase with Q^, at higher the ratio derived 
from HERMES is above unity and they appear to be in- 
dependent of within error bars. 

In Figure 0] we further illustrate the origin of this be- 
havior by plotting separately the numerator (data points) 
and the denominator ( "theoretical" curves representing a 
pQGD based parameterisation), of the ratios Runpoi and 
Rpoi , respectively. We also plot the integral of Fi (dotted 
line) in order to show the effect of both Ai and of the 
Q^-dependent factors that come into play in the defini- 
tion of gi. All quantities are plotted vs. x = (x), i.e. the 



average value of Bjorken x for each spectrum, defined in 
Section II. Notice that the value of (x) increases with Q^. 
The trends in the figure suggest therefore that similarly 
to what observed in DIS, in the resonance region there 
are corrections beyond DGLAP evolution that are posi- 
tive at large x, and negative at smaller x, the threshold 
being defined by: x w 0.33 - 0.43 and ^ 1 GeV^. 
However, while these corrections are comparable in size 
for both the polarised and unpolarised case at large x, at 
low X there seems to be a much larger non-perturbative 
effect for the polariseddata. While possible explanations 
have been suggested e.g. in [23, it is clear that more 
data in this region would help disentangling the de- 
pendence and the possible size of the non-perturbative 
effects. 



B. Comparison with Phenomenological 
Parameter isat ions 

The unpolarised structure function F2 has been eval- 
uated from three different phenomenological fits to DIS 
data 0, I28I |29| and scaled to the same values as 
for and r\''^ to take into account the large effect of 
scaling violation at large x. 

The parameterisation from Ref. |l5l | is a modification 
of the PDF |23| that purports to include target mass ef- 
fects through an effective change of variable, higher twist 
effects at high x, and a factor that enables an extension 
of the fit down to the photoproduction limit. The ALLM 
parameterisation is based on a reggeon and pomeron 
exchanges and it was constructed for the high limit. 
We included it because it can be extended to very low 
values. The NMG parameterisation includes a fit 
of HT terms using world data with W'^ > 10 GeV'^. 

The ratios R^^^^i ^ jrcs/jDis ^.^^ ^ms ^ frcs/pDis^ 

where I^^^ and Vf^^ are the integrals calculated with 
these phenomenological parameterisations, are shown in 
FiglSjin the top and bottom panel, respectively, for sev- 
eral Q-^-values. 

In the unpolarised case the slope of the ratio is less 
evident compared to the previous method shown in Fig|2| 
and the ratio is well consistent with unity for > 2 
GeV^. In the polarised case, the ratio at higher de- 
rived from HERMES data is consistent with unity in- 
side the experimental errors and still independent of Q^. 
However, at low the uncertainty on the parameteri- 
sation is bigger than what found from the PDFs. 

Since these phenomenological parameterisations are 
obtained by fitting deep-inelastic data even in the low 
region, they can implicitely include non-perturbative 
effects and this may explain the "observation of duality" . 
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C. Size of Non-perturbative Contributions 

In order to understand the nature of the remaining 
dependence that cannot be described by NLO pQCD 
evolution, we studied the effect of TMC and LxR on the 
ratios R^^poi ^"^^ ^poT ■ ^he analysis was performed by 
using X as an integration variable, which avoids the am- 
biguities associated to the usage of other ad hoc kinemat- 
ical variables. We used standard input parametrizations 
with initial scale Ql = 1 GeV^. Once TMC and LxR 
have been subtracted from the data, and assuming the 
validity of the twist expansion, Eas. (|3ll9|l in this region, 
one can interpret any remaining discrepancy in terms of 
HTs. 

We notice that although we did not consider NNLO 
calculations, these are not expected to alter substantially 
our extraction since, differently from what seen originally 
in the case of ^3, these have been proven to give a rela- 
tively small contribution to F2 [soj . 

The value of as{M'^) that was used in our calculations 
corresponds to the one given for the DIS parameterisa- 
tions. It has been long noticed that a correlation exists 
between as and the extracted values of the HTs (see [sj 
and references therein and the recent highly accurate de- 
termination in Ref. H^)- It is exactly because of this 
correlation that we keep its value fixed from evaluations 
in a region where the HTs contribution is negligible. This 
statement is equivalent to saying that as cannot be ex- 
tracted reliably from large x data. 

TMC have been evaluated using Eq.l^J for the unpo- 
larised case, and Eas. H4l2U|) for the polarised data. Al- 
though this procedure disregards parton off-shell effects 
that mig ht be important in the resonance region (see 
Refs. p33. l34| ] . we emphasize here its power expansion 
character, and we set as a limiting condition for its va- 
lidity, that the inequality: x^A'P/Q^ < 1, be verified 0. 
Therefore, current treatments of TMC in the resonance 
region are uncertain for values of ^ 1.5 GeV^. 

LxR effects arise formally from terms containing pow- 
ers of ln(l — z), z being the longitudinal variable in the 
evolution equations, that are present in the Wilson coef- 
ficient functions C{z). The latter relate the parton dis- 
tributions to e.g. the structure function F2, according 
to: 



NS 



27r ^ 

9 



dzCNsiz)qNsix/z,( 



(24) 

where we have considered only the non-singlet (NS) con- 
tribution to F2 since only valence quarks distributions 
are relevant in our kinematics. The logarithmic terms in 
Cns{z) become very large at large x, and they need to be 
resummed to all orders in as. This can be accomplished 
by noticing that the correct kinematical variable that de- 
termines the phase space for the radiation of gluons at 
large x, is = Q^{1 - z)/z, instead of 1110. As 
a result, the argument of the strong coupling constant 
becomes z-dependent: as{Q^) as ((3^(1 — z)/z) (see 



|37| and references therein). In this procedure, however, 
an ambiguity is introduced, related to the need of contin- 
uing the value of a s fo r low values of its argument, i.e. 
for z very close to 1 [Sg • The size of this ambiguity could 
be of the same order of the HT corrections. Nevertheless, 
our evaluation is largely free from this problem because 
of the particular kinematical conditions in the resonance 
region. We are in fact studying the structure functions at 
fixed W^, in between 1 < W'^ < 4 GeV^. Consequently 
increases with x. This softens the ambiguity in as, 
and renders our procedure reliable for the extraction of 
HT terms. We illustrate this situation in Fig. 6 where we 
plot the value of as at = 10 GeV^, and we compare 
it with the resummed value in the resonance region, at 
Q'^{l — z)/z for a fixed average Q^, and at Q^(x)(l — z)/z, 
with Q^{x) = W^{x)/{1 - (x)), and (x) = 0.83. 

All of the effects described in this Section are summa- 
rized in the upper panel of Figd In the figure we plot 
the ratio R^npoi^ from Eq.ljHJ, where the numerator is 
obtained from the experimental data, while the denom- 
inator includes the different components of our analy- 
sis, one by one. For unpolarised scattering we find that 
TMC and LxR diminish considerably the space left for 
HT contributions. The contribution of TMC is large at 
the largest values of because these correspond also 
to large x values. Moreover, the effect of TMC is larger 
than the one of LxR. We have excluded from our anal- 
ysis the lowest data point at « 0.4 GeV^ because of 
the high uncertainty in both the pQCD calculation and 
the subtraction of TMC. Also, the pQCD calculations 
at « 1 GeV^ differ from the ones obtained by using 
the available set of parametrizations perhaps because the 
latter are extrapolated well beyond their limit of validity. 

Similarly, in polarised scattering the inclusion of TMC 
and LxR decreases the ratio i?^^ (Fig[7| bottom panel). 
However, in this case these effects are included almost 
completely within the error bars. We conclude that du- 
ality is strongly violated at < 1.7 GeV^. 

The difference between unpolarised and polarised scat- 
tering at low can be attributed e.g. to unmeasured, 
so far, dependent effects, both in the asymmetry, Ai, 
and in 32- Furthermore, a full treatment of the depen- 
dence would require both a more accurate knowledge of 
the ratio R in the resonance region, and a simultaneous 
evaluation of (72- The present mismatch between the un- 
polarised and polarised low behavior might indicate 
that factorization is broken differently for the two pro- 
cesses, and that the universality of partonic descriptions 
no longer holds. 

In Figures |S| El we address explicitely the question of 
the size of the HT corrections. We define them for F2 as: 



H{x,Q^) 
Cht{x) 



Q2 {Fr{x,Q^)-Fi'') 
H{x,Q^) 



7PQCD 



(x/g2) 



Q 



2Fr{x,Q^)-F2 



LT 



F^^ 



(25a) 



(25b) 
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A similar expression is assumed for gi. Cht is the so- 
called factorised form obtained by assuming that the 
dependences of the LT and of the HT parts are similar 
and therefore they cancel out in the ratio. Although the 
anomalous dimensions of the HT part could in principle 
be different, such a discrepancy has not been found so 
far in accurate analyses of DIS data. The HT coefficient, 
Cht has been evaluated for the three cases listed also 
in Fig[7| namely with respect to the NLO pQCD cal- 
culation, to NLO-hTMC and to NLO-hTMC-f LxR. The 
values of 1 + Cht/Q^ are plotted in Fig|Sl(upper panel) 
as a function of the average value of x for each spectrum. 
One can see that the NLO-|-TMC-|-LxR analysis yields 
very small values for Cht in the whole range of x. Fur- 
thermore, the extracted values are consistent with the 
ones obtained in Ref.[^ using a different method, how- 
ever the present extraction method gives more accurate 
results. Because of the increased precision of our analy- 
sis, we are able to disentangle the different effects from 
both TMC and LxR. 

In the polarised case (Fig. 8, lower panel) the HTs are 
small within the given precision, for > 1.7 GeV^, but 
they appear to drop dramatically below zero for lower 
values. The inclusion of TMC and LxR renders these 
terms consistent with zero at the larger values, but 
it does not modify substantially their behavior at lower 

In Figl^ (upper panel) we compare our results in the 
unpolarised case to other current extractions of the same 
quantity. These are: i) the extractions from DIS data, 
performed with the cut: > 10 GeV^ Jsl [ii ^ ; ii) 
the recent DIS evaluation by S. Alekhin ^fll using a cut 
onW^ > 4 GeV^, and including both TMC and NNLO; 
in) the results obtained within a fixed framework in 
Ref. , including both TMC and LxR. We notice that re- 
sults obtained in Ref. ^3 in the deep inelastic region also 
including both TMC and LxR yield small HT coefficients, 
consistent with the ones found in Ref. 0. However, while 
most of the suppression of the HT in the resonance re- 
gion is attributed to TMC, in the contribution of 
TMC is small and the suppression is dominated by LxR. 
In other words, the behavior in the DIS and reso- 
nance regions seems to be dominated by different effects. 
In Fig. 9 (lower panel) we compare the HT coefficients 
in the unpolarised and polarised case. One can notice a 
considerable discrepancy at < 1.7 GeV^. 

Our detailed extraction of both the dependence and 
the HTs in the resonance region establishes a background 
for understanding the transition between partonic and 
hadronic degrees of freedom. In particular, we seem to 
be detecting a region where the twist expansion breaks 
down, and at the same time, the data seem to be still far 
from the limit, where theoretical predictions can 

be made |43j . This breakpoint is marked, for instance, by 
the discrepancy between polarised and unpolarised scat- 
tering at < 1.7 GeV^. More studies addressing this 
region will be pursued in the future, some of which arc 
also mentioned in |4^ . In particular, a breakdown 



of the twist expansion can be interpreted in terms of the 
dominance of multi-parton configurations over single par- 
ton contributions in the scattering process. In order to 
confirm this picture it will be necessary to both extend 
the studies of the twist expansion, including the possible 
dependence of the HT coefficients and terms of order 
0(1/(5^), and to perform duality studies in semi-inclusive 
experiments. Finally, a number of scenarios have been 
studied in [23| , that consider SU (6) quark parton model 
breaking effects and that preserve duality in polarised 
scattering. The results of Ref.|23l are consistent only 
with the larger behavior of the data, although they 
do not address explicitely the question of duality viola- 
tions. 



IV. CONCLUSIONS 

In summary, we presented a study of parton-hadron 
duality in both unpolarised and polarised scattering. The 
latter was obtained by using the first experimental data 
for the polarised structure function of the proton .gf(a;), 
for values larger than 1.7 GeV^. Parton-hadron dual- 
ity was analysed within a QCD context. A pQCD NLO 
analysis including target mass corrections and large x re- 
summation effects was extended to the integrals of both 
unpolarised and polarised structure functions in the res- 
onance region. Within our context, duality is satisfied 
if the pQCD calculations agree with the data, modulo 
higher twist contributions consistent with the twist ex- 
pansion. Although the latter are found to be very small 
for the unpolarised structure function, we do not con- 
clude that parton-hadron duality holds straightforwardly. 
On the contrary, our findings seem to unveil a richer 
dynamics both at a; — > 1 and at small x. This observa- 
tion is substantiated by the fact that duality holds when 
comparing data in the resonance region with phenomeno- 
logical fits which contain some additional "non conven- 
tional" dependence, beyond what predicted by the 
twist expansion. Most importantly, the coefficient of the 
HTs extracted using data in the resonance region only, 
is smaller, and therefore not consistent with the one ex- 
tracted in the DIS region. Finally, while the size of the 
HT contributions is comparable in both polarised and 
unpolarised scattering at larger x and values, at low 
X and we find large negative higher twists only in the 
polarised case. 
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FIG. 1: Ratio between the three different variables x' , ^ and 
xw defined in the text and the Bjorken variable function 
of X. 
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FIG. 2: The integral /„ defined in the text plotted vs. 
(dashed line), compared to the Mellin moments defined in 
Eg. (1101 (full lines). All quantities have been calculated for 
illustration using the parton distributions function parame- 
terisation from Ref. l2a|. 





FIG. 3: Ratio between the integral of the structure function 
as measured in the resonance region and as parameterised in 
the DIS region as a function of . The top panel refer to 
the unpolarised case, while the bottom panel to the polarised 
one. One band represents the theoretical uncertainty due to 
the use of different LO and NLO parameterisations: MRST 
113, CTEQ m, GRV94 [H, GRV98 H. The other band 
represents the experimental uncertainty, that is the sum in 
quadrature of the statistical and systematic uncertainties of 
the data in the resonance region. 
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FIG. 4: The integrals F''" , Eq.lgl, and r^", Eq.© plotted 
vs. the average value of Bjorken x defined in the text, ob- 
tained using the data from Refs. 1^. Holl . Experimental data 
are compared to the integrals of the valence component of 
the structure functions gi (full line), F2 (dashed line), and Fi 
(dotted line), calculated using NLO parameterisations. 
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FIG. 5; Ratio between the integral of the structure function 
measured in the resonance region and as parameterised in the 
DIS region [THHIH as a function of Q^. The top panel 
refers to the unpolarised case, while the bottom panel to the 
polarised one. Notations are as in Fig|3] 
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FIG. 6: as calculated to NLO for different forms of its ar- 
gument: at = 10 GeV^ (dashes), at Q^(l - z)/z 
(full line), and ^ W'^x/{1 - x){l - z)/z (dot-dashed line). 
Calculations including LxR in the resonance region use the 
latter form. 
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FIG. 7: Ratio between the integrals of the measured structure 
functions and tlie calculated ones plotted as a function of Q^. 
The calculation includes one by one the effects of NLO pQCD 
(squares), TMC (open circles) and LxR (triangles), The top 
panel refers to the unpolarised case, while the bottom panel 
to the polarised one. 
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FIG. 8: HT coefBcients extracted in the resonance region ac- 
cording to EQ. II25al . Shown in the figure is the quantity: 
1 + Cht {x)/Q^ . The top panel refers to the unpolarised case, 
where we show the HT term obtained by considering only 
the NLO calculation (squares); the effect subtracting TMC 
(open circles); and the effects of subtracting both TMC and 
LxR (triangles). We show for comparison the values obtained 
from the coefficient H obtained in Ref. lsoll using DIS data and 
including the effect of TMC. The bottom panel refers to the 
HT coefficient in the polarised case. Notations are the same 
as in the upper panel. 
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FIG. 9: Comparison of the HT coefficient displayed in Fig. 8, 
witfi other extractions (upper panel). The triangles and 
squares are the same as in Fig. 8 and they represent our deter- 
mination in the resonance region. Our results are compared 
with extractions using DIS data only. The striped hatched 
area corresponds to the early extraction of Ref . 13911 . The full 
dots are the central values of the extractions in Refs. lioll and 
IZH . These are compared with the more recent extraction 
of Ref. 113 which includes also TMC. Results obtained in the 
resonance region, in the fixed analysis of Ref.0 are also 
shown (stars). In the bottom panel we show the comparison 
between the HTs in the resonance region, in the polarised 
(triangles) and in the unpolarised (stars) cases. 



